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Forces between particles 


Essential to the simulation of stochastic models of particles is the 
consideration of forces acting on each particle (cf Schumacher, 1996). 
Since in theory the forces acting between two particles extend their 
effects to infinity, one has to make approximations. The degree of justi- 
fication to these approximations depends on how the force in question 
vary with distance. 
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Figure 1 Coulomb force vs distance in log-log scale 


In the case of the Coulomb force where the potential energy is 
u(r) = qiqo/(4ae0e,r) and the force f(r) = —du(r)/dr is 


f(r) = u@/(Areverr’), 


the justification is high as is seen in Figure 1 where the force-distance 
is a decreasing straight line on the log-log graph, gq; and q are protons, 
and r the distance in metre. 
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The effects of the Coulomb force are prominent when sizes of the 
particles are small. In quantum mechanics where the scale is atomic, 
for instance, the problem becomes one in which there are a large num- 
ber of particles interacting with one another, i.e. an n-particle problem. 
In theory the Hamiltonian operator can be applied and then the prob- 
lem solved numerically. But in general this is not possible due to the 
too many particles involved in the calculation. Theoretical solution 
is possible by the various methods of approximation, for example the 
Hartree method (cf Brown, 1972) where the best wave function is found 
in terms of the one electron function, i.e. orbitals, ® or a Hatree-Fock 
approximation which reduces the number of equations to n/2. 


Problem definition 


Walls in a 3-d Voronoi tessellation are isomorphic to bonds that 
link between its cells. Therefore we can redefine the problem of a 
particle passing through the hole in a wall, into the cell chamber and 
out through a hole in another wall, and so on, as that of the passage 
of a mathematical particle through a tube that links between two cells, 
into a cell and out through one of the other tubes. In reality the solid 
parts which make up the partitions have thickness, therefore the holes 
through the walls will have a nonzero thickness and the volume of the 
cell chamber is smaller than that of the otherwise mathematical Voronoi 
cell. 


Each vertex of a Voronoi polygon has three walls of the latter at- 
tached to it. The centre of gravity {+ of these three faces is calculated 
and linked together. This truncates the coigns and produces for each 
polygon its dual self. The next step is to link together the midpoints 
of all those edges of the polygon that have a vertex in common. The 
edges of this last polygon bound and define the void volume of the cell 
chamber. 


Next, the cross section of the hole within each wall is taken to be 
the area on the original wall which is bound on all sides by the second 
order covering lattice of that wall. 


The filtering membrane in this case is assumed to be isomorphic 
and homogeneous. The size of the particles is taken to be reasonably 
smaller than the void in general, so that the attraction between parti- 
cles and the adsorption to the wall play a more prominent role than 


} aka centre of mass, centroid. 
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the physical blockages by individual particles. The particles are all 
assumed to be of the same size. 


Algorithm 1 is an algorithm to do filtration that is being developed. 
Here N contains the neighbourhood information 


Algorithm 1 Filtration in Voronoi tessellation 


find (v,,Ca) <_ Voronoi tessellation; 

find d; + Delaunay triangulation; 

find N, < cell neighbours; 

¢ + Cq which lie completely within [0, 1]°; 


Uv € U, which belong to some c; 


My work on filtering membranes results in another Voronoi algo- 
rithm and programme which is different from the programmes for the 
same purpose used elsewhere. Because it is also written anew from 
scratch, this new algorithm has nothing to do with the previous two as 
regarding the data structure and the logic of its method. 


The current convention for variables is this. A single alphabet or 
entity means a list, for example c is the cell list and va is the list of all 
vertices original created. Another example is b, a structure for bonds 
which contains the list of the cells of bonds, the number of vertices 
of the face represented by each bond, the list and the map of these 
vertices, the ordered list of the vertices of each face and another list 
similar to this but cyclic, and the number of cells connected to each 
bond. Similar to the structure of b is that of bdr, or bg in Algorithm 2, 
which is a list of the bonds along the border, each of which is connected 
to only one cell. 


Two entities zy makes x of y, for instance vc is vertices of cells, 
which is a structure that contains the number of vertices of each cell 
and the list and the map of all these vertices. I try to preserve the space 
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by creating a short but easy to understand naming convention. Also, 
the lines are put together, which means that the number of lines of the 
codes is more than what appears here, some of the lines having six 
logical lines or more to them. But the structure of the programme still 
remains intact, therefore it should be possible without much difficulty 
to compare the programme, which shall be published elsewhere, with 
Algorithm 2 which explains it. 


Three entities xy” or y? are yxy matrix of x, where the latter relates 
two entities of the former, for example bcc is the bonds mapped on to 
a cell matrix. In the algorithm, v, and cp, are respectively the vin and 
cin in the programme, the vn and cn there being verbosely described 
in the algorithm as the number of v and c. 


Algorithm 2 Voronoi data structure for the study of membrane filters. 


(vq, v2) + find a 3-d Voronoi tessellation; 

Un + index of 0 < vu, < 1; 

Cn < index of c all the v’s of which are in vy; 

UV Ua E Un} 

Uc — VE E Cnj 

t, + find a 3-d Delaunay tessellation; 

for all t; € t, which are connected to two cells do 
te tj 

else 
ba « ti 

end 

beget 


be ba; 
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find Delaunay triangulation in (3 — 1) dimensions for all faces; 


order the vertices of each face; Oo 


Molecular dynamics 


The Lennard-Jones function can be written u(r) = —A/r° + B/r”, 
where A and B determine respectively the attractive and repulsive 
parts. Let the range parameter o = (B/A)'/° and the energy param- 
eter ¢ = A?/4B, and the equation becomes u(r) = 4¢ [(a/r)!* — (a/r)°]. 
In other words, o is the diameter of one of the atoms and é« is the 
well depth, i.e. energy constant. The force here is not F = —du/dt = 
(24¢/r?) [2(a/r)"? — (a/r)°], but F = —du/dr. 


Another approach in molecular dynamics simulation is to use nu- 
merical methods on the Newton’s equation of motion. The most widely 
used is perhaps the Verlet algorithm, shown here as Algorithm 3. The 
initial values are rp and 11. 


Algorithm 3 Verlet algorithm, cf L. Verlet (1967) 


fori=1tondo 
find fi; 
rig © 2r5 — 74-1 + fit? /m + O(At*) and 
Ui & (Ti41 — Ti_-1)/(2At) + O(At?). 


endfor n 


There are variants and modifications of Algorithm 3, for example 
the Leapfrog Verlet algorithm which has three steps instead of two, 
that is Un41/2 = Un—1/2 + fn/mAt + O(AG), rng = tn + Ungi/2At + O(At*) 
and vp = (Vadis. + Va—1j2)/2 + O(At*); or the velocity Verlet algo- 
rithm where rp4yi = Tn + UnAt + frAt?/(2m) + O(At) and Uns1 = Un + 
At(fnai + fn)/(2m) + O(A#), or in the form normally used in practice 
in which vpiiy2 = Un + frAt/(2m), ta4r = Tn + Ungi2At and unsyi = 
Un + fr+idt/(2m). 
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The van der Waals force is an attractive force acting between mole- 
cules. In the case of gases, each molecule consumes some space, so the 
dynamic volume is less than overall volume by an amount bn when b 
is a constant and n the number of molecules. Here } is Nav, where v is 
the volume slightly larger than the volume of each molecule of gas and 
Na the Avogadro number, N4 = 6.022 x 10-7? mol~!. The pressure of 
gas we see is the pressure of gas detected, which is less than the real 
pressure by F/A, where F = >> Fj and F; = > f;. Here i runs from 1 to 
n and j from 1 to (n—1) in the same set of molecules. The reduction of 
the total force comes from all molecules, but this reduction from each 
molecule is in turn affected by those molecules around it. Therefore 
i # j and j runs from 1 to (n — 1) as mentioned above. Suppose that 
each molecule attracts another molecule by a force a. Then, in a given 
volume V, F = R= 0,0; fig = Ulm — a/V = n(n — 1)a?/V?. For 
very large n we can say that (n—1)  n, and therefore (P+an?/V)(V — 
bn) = nRT. 


The nature of this mutual attraction between molecules, which re- 
duces the pressure in gas, comes to light in the case of solids. The 
distribution of charges around a neutral atom in solid fluctuates in the 
time scale r < 107! s. This charge imbalance makes each atom behave 
as an electric dipole. This electric dipole has an electric field around it, 
which affects other atoms nearby and binds the two together. 


The electric dipole moment is p = ga, where a is a vector from the 
negative to the positive charge. The electric field around a charge q; has 
a magnitude E = q/(47eor”). It acts on a nearby charge q with a force 
of magnitude F = u = mE = mq@/(4neor), where u is the potential 
energy of the two charges. This force itself is the Coulomb force, and 
the electric potential around q is u/q., that is V = q1/(47€o0r). 


For an electric dipole, 


Ve 1 qq @ wet) 

(4n€0) il r2 (4néorir2) 
Since a < r, a being of the order of the atomic size, rir. = r*. Let 
@ be the angle between a and r. Then, also since a < r, we have 
r2 —1 & acos@. Therefore we now have V & qgacos6/(4meor7) = 


pcos@/47eor? and consequently E, = —OV/Or = 2pcos6/(4zeor*) and 
Eg = —OV/(r00) = psin6/(4meor°). 


We neglect Eg for long-range interaction, and simplify E, to E, = 
a/r°. We shall hereafter use the terms dipole and atom interchangeably. 
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If this electric field is produced by atom p;, it could induce in another 
atom p2 = ak£,, where a is the molecular polarisability of the second 
atom. This second atom will have an energy of interaction with the first 
one, u = —poE, = —aE? = —aa’/r°. This force is attractive. 


The repulsive force is more complicated, and is generally thought 
to arise from the Pauli exclusion principle and the Coulombic repulsion 
of the electrons in the outer orbit. It is generally assumed to be either 
u = A/r’? or u = aexp(—r/p), where p is a range parameter (cf de 
Podesta, 2002). 


The Lennard-Jones potential is often written as 


u = —4e [(a/r)° — (a/r)"*] 


where o is a range parameter that indicates the approximate size of 
an atom and ¢€ an energy parameter that indicates the strength of the 
interaction between atoms. In other words, it is the minimum value of 
U. 


Figure 2 Lennard-Jones potential 


The Lennard-Jones potential between two atoms has the minimum 
value u = —€ at r = 1.12250. Figure 2 is a plot y = —4(1/x° — 1/z'%), 
where y = u/e and z = r/o. The force which the potential curve of Fig- 
ure 2 acts on another particle is shown in Figure 3. The force produced 
by a dipole is calculated from F = —du/dr = 24(e/o)[(o/r)’ — 2(¢/r)'9). 
But in our units y = u/e and x = r/o above the equation becomes 
y = 24(1/a’ — 2/2x'), which is plotted in Figure 3. 
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Figure 3 Force corresponding to the Lennard-Jones potential 


Figure 3 shows that the force is zero at the distance r = 1.12250 
away from the particle. Closer than this point the repulsive force in- 
creases rapidly towards infinity. Further away from this point, however, 
the attractive force increases to a maximum and then gradually dies 
down. This maximum force occurs when d?y/dz? = 0, where y = u/e 
and « = r/o, that is at r = 1.24450.: 


In gases, the molecules travel past one another so fast that they 
hardly notice the wells of negative energy surrounding other molecules, 
let alone stop and rest there. But even here the attractive force produced 
by these wells is probably what give rise to the term an”/V? in the van 
der Waals equation. Solid molecules, in contrast, hardly have kinetic 
energy and therefore prefer to sit in such wells of their neighbours, 
which is the reason why solids hold and neither flow as liquid does 
nor disperse like gas. 


The most cohesive structures of solids are crystals, where the 
Lennard-Jones potential culminates in a minimum cohesive energy of 
the lattice. This cohesive energy is U = Nauj/2. Here u; is the sum- 
mation over all pair potential energies of Lennard-Jones type, uj = 
pees u(rj;). Let ri; = aijro, where ro is the nearest neighbour distance. 
Then U = —2eN,[(o/1r0) Ao — (0/70) Ai2], with the lattice sums being 
Ap = Y(1/a§;) and Ary = 3(1/aj7). These lattice sums are calculated 
from the network structure in terms of infinite series that converge very 
quickly. 
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To find the optimum value of a/ro we shall let b = o/ro and write 
differentiate dU/db = —2eN,4(6b° Ao — 12b!'Aj2) = 0. As a result, ro = 
a(2A12/A6)'/°, the cohesive energy per mole U = —(A?2/2Ai2)Nae and 
the cohesive energy per molecule u = U/N4 = —(A2/2Ai2)e. 


The first algorithm we shall now develop is a reasonably simplified 
one. But it turns out to my surprise that this simple algorithm links us 
to the continuum percolation of hard spheres or even of particles with 
irregular shapes. This is because of the assumptions that we shall make. 


One such assumption is that the variance of the filter cell size 
distribution is small. For this, Jafferali (1995) would probably have 
applied his favourite constraint code on a Voronoi tessellation. The 
approach which I developed uses the operator G(-) on a VY. It is more 
logical and therefore is what we shall use for now. 

This centre of gravity operator G(-) has been described (Tiyapan, 
2004) together with the algorithm for the case of two dimensions. In 
the present three-dimensional case for each polyhedron we find c.g.’s 
of all its tetrahedra instead of polygon those of all its triangles. Once 
all these centres of gravity have been found we can then add them all 
up in a similar manner done before. Algorithm 4 tells us how this is 
done. 


Algorithm 4 Centre of gravity of three-dimensional Voronoi Tessellation 


for each polyhedron do 
find the centre, which is the average coordinates of all the vertices; 
for all faces do 
triangulate the face; 
for each triangle do 
find centre of gravity of the tetrahedron made by 
this triangle and the centre; 


endfor 
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endfor 


find the centre of gravity of the polyhedron from the c.g.’s of all 
its tetrahedra; 


endfor _ 


One favourite constraint that has been related to a Voronoi Tessel- 
lation in application is that which limits the minimum distance between 
any pair of the nuclei. The reason given for this is usually that in nature 
cells are of similar size and shape. This implies that the nuclei never 
get together closer than a certain distance, which can be found empiri- 
cally by physically making some observations and measurements. The 
cause of this is often put to the hydrostatic pressure inside the cells, 
which pushing at the cell wall against the pressure outside at the same 
time pushes the nucleus back, away from them. For this reason a rule 
of thumb is created as a criterion used in the generation of Voronoi 
tessellations for the study of various kinds of structure. 


Centre of gravity is a very important point for each object. The 
stability of a double-decker is determined not by its height but by the 
location of the c.g. All forces acting on an object with no angular ve- 
locity can be represented by an imaginary force passing through the 
c.g. In martial arts when it comes to a face-to-face combat you want 
to keep your eyes on the c.g. of your opponent. The limbs and any- 
thing may fly all over the place but wherever his centre of gravity be 
ultimately that is his move. Though in this last case one has to keep in 
mind the body is nonrigid, and each of its components has a different 
momentum. 


In doing the simulation, a Voronoi tessellation is first created using 
Poisson points as nuclei. Then for each cell the centre of gravity is 
found. And then another Voronoi tessellation is created using all these 
centres of gravity calculated as nuclei. Repeating this procedure, the 
centre of gravity function G(-), n times eventually gives us the Voronoi 
tessellation V(G"(x)) we want. Here x is the set of all the nuclei. 


The use of centre of gravity this way is new and has not been found 
in literature. However it is simple while at the same time provides 
a logical starting point, a reasonable basis to work upon which can 
satisfactorily replace the arbitrary limiting distance criteria in use. 
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The second assumption is that the particle size is reasonably small- 
er than the size of the smallest void of the original system. This is the 
assumption which seems to link us to the percolation of spheres in 
continuum, the space considered being that of each void. 


Also, we assume that the particles are attracted towards one an- 
other and towards the walls by various interactions which may include 
the van der Waals force and the force from electrostatic interactions. 


Let @ be the angle that the line from the c.g. of a each cell to the 
mid point of a face makes with the horizontal plane. Then another 
possible assumption is that all particles prefer a trajectory which goes 
through a face which has the maximum 6. 


The third assumption is that the volume of each membrane pore is 
80 per cent the volume of the corresponding Voronoi pore. The fourth 
assumption, blockages due to the clustering caused by the attrition 
forces between particles can occur within pores. This excludes cake 
formation and blockage at entrances to, or exits from the pore. 


The last assumption above is similar to assuming that clustering 
due to attrition can occur in stagnation regions. Since the flow of the 
suspension reaching the membrane is strong, i.e. the pressure high, 
there can be no clustering there. The flow from one pore into another 
is also faster than the flow within one. This is amount to assuming 
that the effective diameter of the hole within each face of a pore is 
considerably smaller than its diameter. 


Upon reaching the membrane, each particle in the troop seeks out 
the bond closest to it, and begins its journey through the membrane. 
Within the membrane, it always follow the bond which has the greatest 
gradient available. Therefore the top bonds should map to the bottom 
ones one to one. But because the possible blockages during the course 
of operation, this may not be so and the mapping is instead one to 
many. 


First we consider the case in two dimension, so volumes becomes 
areas and the volumetric flow rate the distance travelled. In our discrete 
time, if v is the flow velocity, n the number of particles, p, the total 
volume ratio of the particles and r their radius, then we have nar? = 
vAt. Assuming that v = 1 ms! and r = 100 pm, then n = 10°p,/a or 


approximately 318p, particles for each time step. 
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Next we will study the suspension in a square box in order to find 
out p,. Because we shall assume that the particles flowing through 
the pores percolate as though there is no flow, we will consider here 
the suspension that is simply contained within our box without moving 
about. I feel that the assumption that particles could percolate that way 
in pores is justified since the flow is in steady state, protected from the 
turbulence outside by all the solid structures making up the walls of 
the membrane. 


The percolation programme first generates random positions of the 
particles, then moves apart those which are too close together so that 
in the end they only touch each other. Particles which are separated 
by a distance less than 1.4r move towards each other until they touch. 
Lastly, touching particles never separate. 


We know that the total volume ratio of the particles has the upper 
limit of 0.9069, because that is the density of the closest packing of 
circles on a plane. 


Packing circles on the plane becomes densest when the circles are 
arranged as a hexagonal lattice with the packing density of 7/2V3. 
Packing of spheres is similarly at its highest density if the arrangement 
is that of the face-centred cubic lattice with the packing density 7/3V2. 


Filtering problem when physical blockage is prominent 


Assuming particles to be spherical, the size of particles to be con- 
stant, and that this size is compatible with the size of the holes in the 
walls and the voids so that physical blocking is responsible for most of 
the blockages in the membrane. The hole in each wall is taken to be 
the polygon which results from recursively finding a covering polygon 
for the face twice. 


Redefine the problem of particles’ passage through walls and voids 
as that of particles travelling along edges of the dual lattice of the 
Voronoi structure, i.e. the Delaunay triangulation. Each of these edges 
corresponds to a face in the original physical lattice. To each edge is 
thus ascribed the details of the cross section of the hole through that 
face, namely the shape and the dimension, as well as the gradient it 
makes with the horizontal plane. Upon reaching a vertex, the particle 
ball will choose the next path, i.e. bond in the dual lattice, which has the 
maximum gradient to pass through. However, if this bond is too small 
or if it is blocked, the particle with choose from among the remaining 
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paths the one which has the maximum gradient, and so forth. If the 
size of the particle is approximately that of the hole it tries to pass 
through, within two per cent of the latter, say, then the particle will 
blind the passage at that point. The difference between blinding and 
blocking is in the degree of tightness that the particle sits in the hole, 
which reflects in the degree of difficulty to remove it by backflushing. 
This degree is not constant but a function of the relative size between 
the two parties involved, i.e. the particle and the hole. 


Filtering with very small particles 


At this point a new filtering algorithm is considered and investi- 
gated. We introduce an assumption that the particles are all of the same 
size which is very small compared with the size of the pores. So there 
is neither cake formation nor blinding by a single particle. Assume 
that the solid particles suspended in a fluid medium, being dragged 
downwards under their own weight. 


Assume that the van der Waals force plays a significant part and 
that these interactions between particles are governed by the Lennard- 
Jones equation. Assuming that each particle is spherical and acts as a 
dipole with the electric dipole p = 10e, where e is the electric charge on 
a proton, e = 1.602 x 10-! C. Then (cf de Podesta, 1996) 


Pp 
ro FA 1 
Aneor? ( ) 


and the energy of interaction between two particles becomes 


2 
ip 
, = ak? = —*_., 2 
eS Os Cyaan @) 


where a is the molecular polarisability of one particle under the influ- 
ence of the other. If we assume that the repulsive force acts in such a 
way that the repulsive energy is u = c/r'*, where c is a constant, the 
Lennard-Jones potential becomes Equation 3. 


2 
_ ap c 
Byers ee (3) 
Comparing Equation 3 to the form 
A B 
eka ee (4) 
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we have A = ap’/((47)e3) and B =c. Or if we compare it to the form 
__ G6 (9% 12 
uy = —4e [(2)9 — (2)"], (5) 


then we have o = (B/A)° = (16cm7e3/(ap’))'/° and ¢ = (A?/4B) = 
a’p*/(4c(47)4e4). Notice that here o and € are simply parameters, not 
the Stephan-Boltmann constant and the dielectric constant, o is a range 
parameter and approximates the size of an atom while ¢ is an energy 
parameter which shows the strength of the interaction between parti- 
cles. At r =o the value of u, is zero, whereas u, has the minimum 
value of —e. 


The force between two particles is Equation 6. 


2: 
paar 56 6a hee (6) 


dr (4m)2e2r7 * 8 


At the equilibrium separation, ro, F = 0 and therefore Equation 6 
yields the minimum distance in Equation 7, 


1 
32cene} \ * 
= {| ——- | .. i: 
ee ( ap? ) 7) 


The equation for the minimum energy is obtained by substituting ro 
from Equation 7 into Equation 3, which gives 


2 
2 2 
= - ate. (11) 


Returning to Equation 6, the attractive force has the maximum 
value at the point where dF/dr = 0. Differentiating F in Equation 6 
with respect to r, we arrive at Equation 12. 


dF 42ap* 156c 


es 12 
dr = (4m)efr8 or 2) 


Vaen Sryayudhya, Editor April 2006 171 


Tyabandha Journal of Arts and Science Vol. 3, No. 2 


From this the distance where this maximum force occurs is given by 
Equation 13, 
eas 156c(47)e4 = Al6cne% 
oe Aap? Jap — 
Then by putting Equation 13 into Equation 6 we have this maximum 
force in Equation 16. 


(13) 


6ap* 12¢ 
Fs Oe a ee 14 
(4m)Pegr7 + 7a m) 
ee ( —6ap? —-12c(42) ap? ) (15) 
~ 7 \ (4n)2e2  156e(4r)2e7 
Sap? 
= -—— 16 
52r7 02h (16) 


v(0)=0 


a 


(b) mg 


Figure 4 Force balance with one particle fixed 


Let each of the spherical particles has a radius r = 1 ym and the 
density p = 3,000 kg- m3 (compare the density of silicon, which is 
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2,329 kg-m~%). Then the mass of the particle is m = pV = 3000 x 
(4/3) x (10-6)? = 1.26 x 10-14 kg and the weight is mg = 1.26 x 10-4 x 
9.8 = 1.23x 10-8 N. We may now find the attractive distance, the 
maximum distance whereby two particles will come together under the 
van der Waals force. Figure 4 shows the capturing of one particle by 
another when one particle is fixed in space and another particle has zero 
velocity but is free to move. Recall that the capture occurs when the 
acceleration due to the weight of the particle g equals the acceleration 
a due to f. 


If we suppose that our particle have the same polarisability as that 
of a benzene in its gaseous state, CoH¢, i.e. a = 11.61 x 10-49, then in 
this case a = 11.61 x 10-*° F-!m?# (cf de Podesta, 1996) and Equation 16 
gives us the separating distance between the two particles which gives 
the maximum attractive van der Waals force, which is 444 wm. Giving 
our particle other values of a, with the value of a for methanol gas 
CH30OH, i.e. a = 3.860 x 10-49 F-'m4, we have this distance r = 379 ym, 
and with a = 1.647 x 10-*° F-!m’, that of water vapour, the distance 
becomes r = 336 pm. 


These values are rather large compared with the radius of the par- 
ticles, therefore we shall opt instead to a bigger size of particles, when 
ry = 5 pm. With the radius of five microns, the particle has a mass of 
1.571 x 10-'* kg and its weight becomes 1.64 10°!’ N. Then if we 
adopt the value of a for water vapour, a = 1.647 x 10-40 F-!m4, then 
from Equation 16 the distance where the force is maximum becomes 
r = 167 pm. 


But this is not everything. So far we have only considered what 
two particles will do when one of them is fixed and the other one has no 
initial velocity. There are two other things that can happen, the particles 
may both be moving down beside each other under gravity and the 
path of the particle which is effected by their mutual attraction. Even 
when two particles do not come together, their path can be deviated 
by the van der Waals force. But here for simplicity we shall neglect 
this and assume that particles either come and stay together or they 
experience no mutual force whatever. Only if they come together will 
their final path and velocity be affected, and these depend on their 
combined momentum. 


When particles move relative to each other, their captive velocity 
depends on their relative velocity. They tend to join each other more 
easily if their paths are along side each other and goes in the same 
direction. In the extreme case where both particles move in the same 
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direction with zero relative velocity, they would come together across 
a vast gap in between indeed, had there been no frictional force due to 
viscosity that acts to drag them. 


In the present study, the friction due to the fluid in the medium 
is neglected, together with the relative velocity between the particles, 
for the purpose of deciding whether particles will come together. It 
is expected that on average particles coming to within the capturing 
radius of each other will come together. This radius is larger than the 
radius of maximum force in Equation 16 above. In other words, we 
expect our particles to behave like solid spheres with a well defined 
boundary for their sphere of influence and a much simplified capturing 
mechanism. 


We shall define this radius of captivity, r-, to be one third that of 
the radius of maximum force r,,. From Figure 2 this is the point where 
re = 1.50, and because r,, is here 1.2445 we have r./rm * 1.2. For 
Figure 2, this is the same as saying that r./ro = 1.34. 


One additional point is that, in real situation when particles form a 
cluster their collective value of the molecular polarisability will change, 
and this value will not be the same for the cluster as it is for the in- 
dividual particles. But here we assume that they are the same, which 
means that the capturing radius of a cluster will be the same as that of 
the particles which form it. 


It may worth mentioning here that the molecular polarisability is 
related to an optical property, viz. the refractive index n,, by the rela- 
tionship a = €o(n2 —1)/n, where n is the number density of molecules. 
For gases n = P/(kgT), whereas for liquids n = Nap/m, m being the 
mass. 


Next we shall concern ourselves with clusters of particles thus 
formed. All clusters will be assumed to be a closest-conglomerate of 
particles which form them, with the densest packing density possible 
in three dimensions, that is 7/32 of the face-centred cubic lattice. 
Numerically this is 0.7405. 


Notice that the packing density in two dimensions can be higher 
than this. Packing circles on a plane is the densest of all with its packing 
density of 7/2V3. 
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(a) 


(b) 


Figure 5 The packing density calculation of the closest-packed densities 


Shown in Figure 5 are the hexagonal lattice packing in two dimen- 
sions and the cubic close pack. In the cubic close packing spheres in 
every third layer lie vertically straight on top of one another. Each face 
of its cubic section looks like the second picture in Figure 5. Similar 
to the cubic close packing is the hexagonal close packing spheres in 
every alternate layer of which lie over one another. Both the cubic- 
and the hexagonal close packing have the same packing density which, 
in the case of the former, is calculated, from the second picture, as 
p = >d5us/>>vc, where vs are the total volume of sphere segments 
in the unit cell which has the volume v,. Here ve = (2V2r)> and 
Yous = (8(1/8) + 6(1/2))4ar3/3. In 2 dimensions, the densest packing is 
calculated from the first picture.: 


The density of the packing is }>a,/a,, where a, are the areas of 
the circles and a, is the area of a rectangle. These two areas are namely 
Qc = 2(ar?) and a, = 2r(2V3r), which give the density p = 7/2V3. 
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If the particle size is very small there will be an effect from quan- 
tum mechanics. The de Broglie wavelength, A = h/p, exists for all 
particles large and small. Here p = mv, m = pV and V = rd?/6 for 
a spherical particle. Figure 6 shows a graph between diameter and 
wavelength of particles. 


d/m 


Figure 6 Quantum-mechanical effect in filtration 


We can see that the quantum-mechanical effect when the pore size 
is in the order of micrometre will become prominent if the size of par- 
ticles challenging it is in the order of nanometre. When this is the case 
there will be not only the interference effect of the particle-wave pass- 
ing through slits but also the effect of confined particles which produces 
standing waves of these wave-particles in each pore. These standing 
waves of the particle-wave have A, = 2D/n. 
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